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Abstract 

We show that the intermittent and self-similar fluctuations dis- 
played by a slow crack during the propagation in a heterogeneous 
medium can be quantitatively described by an extension of a classical 
statistical model for fracture. The model yields the correct dynam- 
ical and morphological scaling, and allows to demonstrate that the 
scale invariance originates from the presence of a non-equilibrium, re- 
versible, critical transition which in the presence of dissipation gives 
rise to self organized critical behaviour. 

1 Introduction 

The failure of materials is a complex and complicated process exhibiting 
broad phenomenology. The fracture of heterogeneous media under slow ex- 
ternal loading displays intermittent dynamics and scale invariance, features 
observed in different phenomena involving a huge range of length and time 
scales, from nano-plasticity [20] and micro-fractures [TIJ [36j [231, EH] to earth- 
quakes [l5l [15j |29] . In the propagation of a slow crack, the front advances 
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through a movement that is statistically stationary but characterized by sud- 
den and intermittent self-similar bursts, or avalanches. Bursts area and du- 
ration are power law distributed, and the crack profile is self-affine (see e.g. 
[5j [6] and refs. therein). Besides having practical relevance, this phenomenol- 
ogy, known as "crackling noise" [I2J , characterizes also the dynamics of other 
interfacial phenomena in disordered media, like imbibition, wetting, friction, 
and hysteresis in ferromagnets. 

In analogy with equilibrium phenomena, the absence of typical scales is 
considered a mark of some underlying critical transition [19]. This conjec- 
ture has led in the last decades to the formulation of simplified statistical 
models for fracture [TTJ [27J El E7J with the aim of understanding the origin 
of scale invariance, and identifying possible universalities in different sys- 
tems. However, it is diffuse opinion [?J [38] that these models have until now 
failed in catching what are thought to be the right power laws characterizing 
both the failure of various materials and the seismic behaviour accompany- 
ing earthquakes. Extant models yield exponents significantly higher than 
the experimental observations [TJ, and make the right values "a target for 
theoretical models" [2J. Moreover, they fail in reproducing the microscopic 
statistics in quasi-stationary situations |2j|, making the link with a reversible 
transition problematic. 

On general grounds, observation of scale invariance in many natural phe- 
nomena has led to the concept of Self Organized Criticality (SOC) [31 [T3] . 
in which a system spontaneously sets close to a critical point, in contrast 
with ordinary critical phenomena where some parameter needs to be finely 
adjusted. Attempts have been made to relate the self-similar fluctuations 
observed in fracture phenomena with SOC dynamics [361 I2H HD], and cel- 
lular automata have been devised with the aim of reproducing the power 
laws observed in fractures and earthquakes [2H1 d3l EE2]- SOC models have 
been very helpful and intuitive in illustrating the general properties of the 
critical behaviour; however their matching with real systems is still difficult 
in relation to several points, in particular the correspondence between model 
parameters and real ones |46j. So the critical transition beneath the frac- 
ture phenomenology has remained elusive, even if more evidences have been 
recently emphasized [7J. 

In this paper we show for the first time that the critical behaviour charac- 
terizing the slow fracture and crack propagation in a heterogeneous medium 
originates indeed from the existence of a non-equilibrium critical phase tran- 
sition, separating two distinct phases: an active phase, in which the fracture 
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propagates indefinitely, and a dormant one, in which the system is quiescent. 
We shall introduce an extension of a widely employed statistical model, de- 
scribed in the next section, and use it to reproduce the propagation of a 
planar crack. It will be seen in the following sections that this model yields 
intermittent and self-similar dynamics in quantitative agreement with the 
dynamical and morphological scaling measured in recent experiments [331 E] 
and in simulations [SI ED]- Moreover, a same universal scaling is found to 
characterize the avalanche area, stress, and energy. The nature of the crit- 
ical transition will be discussed in a separate section. Finally we will show 
that dissipation is crucial in generating the observed intermittency by setting 
the systems at the edge of criticality, and connecting the process to the SOC 
phenomenology. 

2 The model 

The model we adopt is an extension of the classical Fiber Bundle Model 
(FBM), which in its original formulation [18J consists of a set of N parallel 
fibers, each one having a random, quenched, breaking strength drawn from 
some identical probability distribution. Applying a slowly increasing load 
parallel to the fibers axis, more and more fibers break and redistribute their 
stress among those still intact, which in turn can break, generating avalanches 
in a domino effect. The probability distribution for the number of fibers 
broken in an avalanche is a power law, with a cutoff at large values. The 
dynamics is non stationary, as the number of intact fibers decreases with 
increasing load, and at a critical stress the system undergoes global and 
irreversible failure. While in its original formulation FBM is a mean field 
model [US [22], several variants have been also devised to take into account 
dimensionality and different stress redistribution rules among fibers [21 [371 
[251 S3 HZl [211 CO- 
In order to describe the propagation of a crack front during its quasi- 
stationary regime, we have complemented the classical FBM with two rules: 
1) Fiber regeneration: Each time a fiber beaks, it is replaced by another one, 
with new random breaking strength and zero initial stress. 2) Energy dissipa- 
tion: When a fiber breaks, part of its elastic energy is dissipated. Thus only 
a fraction of its stress is transferred to other fibers, while a complementary 
fraction is lost. 

These rules characterize the Dissipative Regenerating Fiber Bundle Model 
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Figure 1: Avalanche energy (lower curve) and instantaneous bundle stress 
(upper curve) in the DRFBM as functions of the applied external strain, 
from numerical simulations with N = 10 4 and 5 = 10 -3 . Inset: a snapshot 
of the crack front obtained by the 1-d DRFBM, as described in the section 
discussing the roughness 



(DRFBM) investigated in the present paper, and are motivated by the follow- 
ing considerations: 1 ) Regeneration. The replacement of broken with intact 
fibers mimics the crack propagation, in which regions previously hedged from 
the stress get involved as the overloaded zones yield. The propagation of a 
planar crack, orthogonal to the fibers, can thus be described by a Id version 
of the model. During each avalanche broken fibers are replaced by new ones, 
each failure event representing an onwards local motion of the crack, as will 
be explained in detail in the section dedicated to the roughness. 

2) Dissipation. The dissipation of elastic energy in real fractures may 
happen by several mechanisms: in the Griffith's description of a perfectly 
brittle medium it is due to the cost of opening the crack. In the general case 
other processes can contribute, such as the emission of elastic waves and the 
plastic deformation. Whereas the introduction of dissipation in the standard 
FBM would just cause a slow down of the dynamics, and a delay of the final 
breakdown, we will show that in the DRFBM it is crucial for the onset of 
intermittence in the motion of the crack front. 

We consider here a bundle of harmonic fibers with identical unitary elas- 
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tic modulus, each one having a random, quenched, breaking threshold t 
extracted from an identical probability density p(t). The dynamics of the 
DRFBM can be summarized as follows: starting from zero initial stress, all 
the fibers are subjected to a same slow increase of strain u, until the weakest 
fiber breaks and it is replaced by an intact fiber having zero stress and a 
new random threshold. Soon after, a fraction (1 — 5) of the released stress is 
redistributed within the bundle according to some rule, while the remaining 
fraction 5 is lost. The redistributed stress may cause the breaking of one 
or more other fibers, which in turn will redistribute part of their stress, and 
so on. Different fibers broken from a same redistribution process are all re- 
generated at the same time. In the model we assume a separation between 
the timescales of the internal avalanches dynamics and the driving field rate, 
keeping the external strain fixed during avalanches. An avalanche stops when 
the stress transferred from the broken fibers causes no further failures. At 
this point the bundle is subjected to a new strain increase until a fiber breaks 
and so on, and the dynamics evolves under such an adiabatic driving. It is 
worth noting that each fiber bears a history dependent stress, resulting from 
summing, from the epoch of its generation, the internally redistributed stress 
and the externally applied strain. 



3 The propagation of a planar crack 

Here we consider the DRFBM model in d = 1, with the fibers placed on a 
line and subjected to periodic boundary conditions. Each site is identified by 
its discrete coordinate Xi, to which it is associated a set of successive fibers 
i a . Each time a fiber i a breaks it is replaced by a new one, i a +i, and the 
crack front moves onwards. According to the rules that will be spedified in 
the section devoted to the roughness, this gives rise to an irregular crack 
front profile, like the one sketched in the inset of fig. 1. Since only one fiber 
at time is present at a given site, we shall from now on identify the current 
intact fiber with only the site index i and the related stress with Sj, unless 
otherwise specified. 

In a crack front each point generates stress on other points proportionally 
to the inverse square of the relative distance [22j H3] . Therefore we adopt a 
distance dependent stress redistribution rule |28j such that if the fiber at x% 
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breaks under the stress Sj, the stress of the fiber at Xj increases by an amount 



ASj = K 



for j ^ i and As» = 0, 



(1) 



where X 2 0C j IS taken modulo N, and k is a 5 dependent normalization factor 
assuring that a fraction 5 of is dissipated: ^ Asj = (1 — 5)§i. 

Figure 1 shows the behaviour of the 1-d DRFBM at increasing applied 
strain, obtained from the simulation of a system with N = 10 4 , 5 = 10~ 3 , 
and random fiber thresholds t extracted uniformly in [0, 1). The lower curve 
represents the energy of avalanches E = where the sum is over all 

the failures in the avalanche. We observe that at low strain there are only 
a few rare and small failure events, but eventually a large avalanche starts, 
leading the systems to a state where avalanches of any size occur. At that 
stage each increase in stress due to the externally applied strain balances in 
average the dissipated stress, and the total bundle stress a = Y^iSi/N reaches 
a statistically stationary value < a >, as shown by the upper curve in the 
figure. 

We have extensively investigated the 1-d DRFBM with different values 
of dissipation 5 and size N. For 5^0, under increasing strain the system 
reaches a state in which the bundle stress a is statistically stationary. The 
case 5 = is special and will be discussed in a devoted section. Besides 
the energy E, we have computed several quantities for each avalanche, and 
their probability distributions : The area A, defined as the total number 
of failed fibers; the duration T, as the number of regeneration processes 
occurring within the avalanche, counting as sigle process the simultaneous 
regeneration of more fibers; the stress S, as the sum of the stress born by 
all the fibers broken in the avalanche. In addition, since isolated avalanche 
clusters can take place at the same time in different regions of the crack front 
[3 E21 ED], we have also considered the statistics of the area C of isolated 
clusters, i.e. of locally connected but mutually disconnected failures in the 
same avalanche. 

It is seen that the bursts display self-similar features in all the quantities 
considered. For large enough avalanches the distribution for each quantity, 
y, is well described by 



The exponents r y do not depend on the dissipation rate 5, which instead 
determines the cut-off values y a . By assuming a power law dependence, y Q ~ 




(2) 
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Table 1: Dynamical scaling exponents for a planar crack from the 1-d 
DRFBM 

T if) (j) 

C 1.5 1.4 0.7 
A 1.2 1.3 1.0 
T 1.5 0.7 0.4 



5~^ v , the distributions of a same quantity for different 5 can be collapsed onto 
a unique curve. Figure 2 shows, on double logarithmic scale, the collapsed 
probability distributions of (a) the cluster avalanche area C, (b) the total 
avalanche area A, (c) the avalanche duration T, as obtained from numerical 
simulations of the 1-d DRFBM with N = 10 4 and for different values of 5. 

By means of the data collapse we have got as best estimates for the 
exponents: ta = 12 and tq = 1.5 for the total and cluster area, respectively, 
and tt = 1.5 for the duration. These values, resumed in the table I, strictly 
agree with those found in the experiments [5j [33], in linear fracture models 
[7J with following rediscussions [5j El [30] , and in depinning models (see e.g. 
[7J and refs. therein). The rescaling cut-off exponents ip y are also reported 
in the table I. Since it is expected that the average avalanche obeys the 
scaling < y >m 5^^ y , the exponents must be consistent with the relation 
4> y = 4> y (2 — Ty), that is indeed the case. 

We have also computed the distribution of the energy E and stress S, for 
both total and cluster avalanches, and we have found that they are described 
by the same exponents of the respective avalanche areas tq and ta, indicating 
that the scaling of the size fluctuations is a universal feature, irrespective of 
the measured quantity. 

Finally we have checked that all results do not depend on the probabil- 
ity distribution p(t) from which the fiber threshold strenghts are extracted, 
provided that this vanishes rapidly enough as t — > oo. Finite size effects are 
not generally displayed as they appear only for very small values of N ■ 8. 

4 The roughness of the front 

The self-similar features of crack profiles have been subject of study since 
when they were first observed in the surface of fractured metals [34]. A wide 
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Figure 2: Collapsed probability distributions for the avalanches in the planar 
crack front dynamics as resulting from the DRFBM, relative to: (a) cluster 
avalanche area C; (b) total area A; (c) duration T. Relevant exponents are 
reported in table 1. 



amount of work has then been devoted to characterizing different morpho- 
logical scalings, clarifying their origin and looking for classes of universality 
P,IH]. It also required time to explain some apparent contradictions between 
different experiments, and discrepancies bewtween some experimental results 
and models for elastic depinning (see e.g. [El Ej), but it seems now rather 
well established that different scaling regimes can be attributed to the in- 
terplay between elasticity, disorder and possible spatial correlations in the 
disorder strength [31], HQ] . 

In order to define the position of the crack front in the 1-d DRFBM, we 
assume that each time a fiber yields, the crack advances locally - perpendic- 
ularly to the fibers and to the distance proportional to the stress s 
stored in the broken fiber. Since fibers are harmonic, the total perpendicular 
displacement at Xi is given by = Yl a s ai where a, counts all the consecutive 
failures occurred at site i up to the current epoch. 

The scaling properties of the roughness can be evaluated through different 
quantities. Here we compute 

R] = ((z l+e - Zi f) (3) 
that is expected to scale as £ 2 ^. For the 1-d DRFBM we find a non trivial 
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Figure 3: Scaling of the crack front roughness R in the 1-d DRFBM. Curves 
from systems with different dissipation have been collapsed by assuming the 
same dependence of the upper cut-off on 5. The dashed line has slope ( = 
0.35. 



power law behaviour of Rg. Figure 3 shows the values of Ri as function of 
I for uniform initial distribution of the fiber strength, iV = 2 16 and different 
dissipations 5. The crack profile is self-similar in a wide range and it can 
be noticed that even in this case the dissipation plays the fundamental role 
of setting the upper cut-off £ Q of the scaling range. Different curves corre- 
sponding to different values of dissipation can in fact be collapsed assuming 
the scaling ([2]) with a dependence of the form l a ~ 5~^ 1 . A best fit yields 
( = 0.35 in agreement with those found in experiments and simulations for 
uncorrelated and not too strong disorder [^Ql EI]. We also find ipg = 1.00. 

Very recently, the transition betweeen two different scaling regimes has 
been also observed in a different FBM based model [?] in which the fibers 
interact through a 2-d kernel. The model shows that, as the Young modu- 
lus increases, the front becomes rougher and overhangs and damage islands 
appear, the exponent in the soft elasticity regime being in agreement with 
that observed in the DRFBM. 
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5 The critical transition and the effect of the 
dissipation 

As discussed in the introduction, a major point is to establish whether a 
critical transition is at the origin of the scale free fluctuations observed in 
the propagation of the crack front. In fact, linking irreversible intermittent 
phenomena to a reversible transition may be of help in the identification of 
different universality classes and of key working mechanisms. 

In a phase transition, a system must exist in at least two different states 
separated by a critical point. In order to address this issue we have investi- 
gated the system behaviour with null dissipation (5 = 0). Note that in this 
case all the energy fed remains in the system, and the bundle stress a = J2i s i 
equals the totally applied strain. 

As usual, an external strain was supplied in a quasi-static mode, starting 
with zero initial stress. The 5 = case displayed initially a transient phase 
similar to the case 5^0. As the bundle stress was increased, larger and 
larger avalanches appeared, until eventually the bundle stress a reached a 
maximum (like in fig.l for 5^0) and a very large avalanche was triggered. 
While in the case 5^0 this peak is necessarily followed by a decay and a 
sets at a lower stationary value, in the 5 = case the domino effect provokes 
a neverending avalanche, and the system remains indefinitely in the active 
state, the redistributed stress bouncing among the fibers. Hence, for null dis- 
sipation, the system is driven to a critical point where the crack propagation 
sustains indefinitely, whilst in the presence of dissipation, however small, it 
goes back below the critical point, and only new a stress increase will trigger 
more avalanches. 

We have simulated the 1-d DRFBM with no dissipation and different 
sizes, and computed the fraction of active fibers p, i.e. the average number 
of fibers simultaneously breaking during avalanches, at varying bundle stress. 
The results are shown in fig. 4, where it is seen that close to a critical value 
a c the activity suddenly departs from zero and monotonically increases. We 
have also computed the variance of p, that is proportional to the system 
susceptibility. It is seen to display a peak close to a c , and to increase linearly 
with the system size, as expected from finite size scaling in d=l. 

These features are the hallmark of a critical transition and show that for 
5^0 the system is at the lower edge of criticality. In fact, in this case, the 
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peak stress attained before the stationary state, visible in fig. 1., corresponds 
to the critical value a c . 

Finite size scaling study and detailed analysis of this transition will be 
reported elsewhere. We stress here that the resulting picture is the one 
described by Dickman et al. [19] for the onset of Self Organized Criticality 
(SOC), where the system oscillates continuously between an active phase, to 
which it is driven by the external field, and an absorbing phase where it falls 
because of the dissipation. This demonstates in a clear and definite way that 
crack propagation is a realization of SOC, as also recently conjectured [7]. 

Dissipation, that has been treated as an irrelevant parameter in linear 
elastic fracture models [SJ [TJ [50] and useful regolarizing factor in depinning 
models [32j, turns out to be the key factor for the onset of intermittence, 
establishing a clear link between fracture processes and SOC. In the DRFBM 
it occurs in the bulk, at variance with usual Self-Organized-Critical (SOC) 
models {e.g. sand piles) where it is limited to the boundaries in order to 
obtain the cut-off to disappear when taking the thermodynamic limit. As a 
matter of fact, internal dissipation is hardly unavoidable in real macroscopic 
systems, and there is growing evidence that it governs real manifestations of 
SOC [7J HU HI [35]. It does not affect the exponents and may well be at the 
origin of the large scale cut-off always present in real systems, in alternative, 
or in addition, to finite size effects. Finally, it is also noticeable that, to 
our knowledge, a diverging susceptibility had never been observed in SOC 
models, whereas it characterizes natural phenomena [55] . 

6 Discussion and summary 

We have described the propagation of a slow crack front in a heterogeneous 
materials by means of a stationary Fiber Bundle Model which includes inter- 
nal dissipation and fiber regeneration (DRFBM). In the absence of dissipation 
the system can be in either of two distinct phases, one active and the other 
quiescent, separated by a critical value of the system stress at which the 
susceptibility diverges. The presence of dissipation, however small, drives an 
active system to a slightly subcritical state whereas, by supplying stress, the 
system can be driven from the quiescent to the active phase. The alternance 
of these two processes leads the system to a statistically stationary state 
where the dynamics is scale invariant and intermittent. 

The 1-d DRFBM with the stress redistribution rule eq. (1) yields the 
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Figure 4: Evidences of the critical phase transition displayed by the 1-d 
DRFBM in the absence of dissipation. The fraction of active fibers p suddenly 
departs from zero above a critical value of the bundle stress a. The variance 
a 2 p shows a peak that increases with the system size. 



correct experimental scaling of both dynamical and morphological fluctua- 
tions of the planar crack. Universal statistics characterizes the burst area, 
energy and stress, whereas dissipation sets the upper cutoff to the self-similar 
range. 

The DRFBM establishes a link between SOC models and real systems 
and can be easily extended to describe systems in higher dimension, like 
non planar cracks and earthquakes, or can be employed for implementing 
preexisiting models [H]. It also presents interesting analogies with other 
breaking-healing contact models, recently introduced to describe tribological 
experiments, earthquakes [9] and biological systems [2TJ, that would deserve 
to be investigated. 
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